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ABSTRACT 

Context. A number of pulsating stars with rotational splittings have been observed thanks to the CoRoT and Kepler missions. This 
is particularly true of evolved (sub-giant and giant) stars, and has led various groups to investigate their rotation profiles via different 
methods. 

Aims. We would like to set up some criteria which will help us to know whether a decreasing rotation profile, or one which satisfies 
Rayleigh’s stability criterion, is compatible with a set of observed rotational splittings for a given reference model. 

Methods. We derive inequalities on the rotational splittings using a reformulated version of the equation which relates the splittings 
to the rotation profile and kernels. 

Results. These inequalities are tested out on some simple examples. The first examples show how they are able to reveal when a 
rotation profile is increasing somewhere or inconsistent with Rayleigh’s criterion in a main sequence star, depending on the profile 
and the ( values of the splittings. The next example illustrates how a slight mismatch between an observed evolved star and a reference 
model can lead to erroneous conclusions about the rotation profile. We also show how frequency differences between the star and the 
model, which should normally reveal this mismatch, can be masked by frequency corrections for near-surface effects. 

Key words, stars: oscillations (including pulsations) - stars: rotation - stars: interiors 


1. Introduction 


The CoRoT and Kepler spac e missions have obtained exquisite 
pulsa tion data for many stars (Ba ulin et all 120091. iBorucki et all 
2009). This has enabled the detection of rotational split¬ 


tings in a number of stars, including subgiants and giants 
(Beck et ak L 20 2l Deheuvels et all 1201 2l iMosser et all [2 012. 
Deheuvels et al, 12014) . and a main sequence star dKurtz et al 


20141) . A number of prior studies have also used ground-based 


data to extract rotational splittings. Based on these splittings, the 
above authors have inverted or constrained the differential ro¬ 
tation profile, and hence constrained angular momentum trans¬ 
port within stars. In particular, the core of red giants rotate much 
more slowly than what is expected based on theory, thereby 
pointing to unknown powerful angular transport mechanisms 
whic h operate throug hout the stellar lifetime (Eggenbergeret al., 

l20lllMaraues et alll20HlCeillier et al.Ll2013l). 


There are several inverse methods used to probe the inter¬ 
nal rotation profile of a star. One of these, the regularised least- 
squares (RLS) method, searches for the optimal profile which re¬ 
produces the observed splittings. When carrying out such an in¬ 
version, it is necessary to introduce a priori information. Indeed, 
rotational splittings represent a finite number of measurements 
or constraints on the rotation profile, i. e ., a function defined 
over the interval [0,/?], where R is the stellar radius. As the 
method’s name suggests, this is typically done through a reg- 
ularisation term which reduces the second order derivative of 
the resultant profile. Even then, the obtained solution is not al- 
wa ys sa tisfactory. Indeed, as can be observed in, e.g., Fig. 15 of 
[Deheuvels et al.l (120121) . the resultant profile may change signs. 
Physically, this would correspond to a star with an “onion” 


type structure with counter-rotating shell(s), as pointed out in 
iDeheuvels et ali (l2014t) . A second, potentially problematic, situ¬ 
ation is when the gradient of the rotation rate becomes positive 
(see, e.g.. Fig. 6 of ICorsico et all [2011 ). Although such a situ¬ 
ation can occur and has been observed both at low latitudes in 
the sun (iScho u et al. il 1998b and in a terminal age main sequence 
A star (IKurtz et all 120141) . it seems unlikely in many cases, es¬ 
pecially in evolved stars which are undergoing core contraction 
and envelope expansion. Therefore, it is important to search for 
inversion methods which impose a positive rotation profile and, 
optionally, one which decreases outwards. Before developing 
such a method, however, it is useful to check beforehand whether 
such a profile is compatible with the observations for the chosen 
reference model. 


In the present paper, we investigate under what conditions it 
is possible to obtain a decreasing rotation profile, or one which 
satisfies Rayleigh’s stability criterion, for a set of observed ro¬ 
tational splittings, and a given reference model. In Section [2] 
we show how these assumptions on the rotation profile lead to 
inequalities on the rotational splittings. Section [3] then contains 
two test cases. The first one shows how such inequalities can 
detect when a rotation profile does not obey Rayleigh’s stabil¬ 
ity criterion, and the second illustrates how these can be used 
to reveal a mismatch between the reference model and the star, 
provided one assumes a decreasing rotation profile. A short dis¬ 
cussion concludes the paper. 
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2. Inequalities on rotational splittings 

In spherically symmetric non-rotating stars, the pulsation modes 
are described by three quantum number: the radial order, n, 
which corresponds to the number of nodes in the radial direc¬ 
tion, the harmonic degree, (, which is the total number of nodal 
lines on the surface, and the azimuthal order, m, which gives 
the number of nodes around the equator. For each pair ( 77 , 7 ) 
there are 27 + 1 modes with the azimuthal order ranging from 
-7 to 7. These modes are degenerate, i.e., they all have the same 
frequency. If, however, the star is rotating, these modes will no 
longer be degenerate. Furthermore, if the rotation profile is slow 
and only depends on the radial coordinate, r , then, based on first 
order effects, these modes will be evenly spaced by a quantity 
known as the rotational splitting, s n /. Using the variational prin¬ 
ciple, it is possible to derive a relation between the rotational 
splitting and the rotation profile, Q (e.g. lAerts et all 2010 ): 

Sn,e = — ~ ^ = (1 - Cn,t) r K n/ (r)Q(r)dr, (1) 
m Jo 

where a> is the pulsation frequency, C„j the Ledoux constant, 
K n { the rotation kernel, and R the stellar radius. In the above 
formula, we have assumed that the star is viewed from an in¬ 
ertial frame (as opposed to a co-rotating frame), and are mak¬ 
ing use of what could be called the “prograde convention”, i.e., 
modes with positive azimuthal orders are prograde. If the op¬ 
posite convention is used, then the splitting is defined as s n j = 
(rVn,l,—m (j hiJ.0 ) I^Tl. 

The Ledoux constant takes on the following expression 
(Ledouxj|1951b : 


Tbfl+'lVdr 

L n ,e = —r -, U) 

f 0 [J 2 + 7(7 + l)? 7 2 ]pr 2 dr 

where £ and r\ are the radial and horizontal Lagrangian dis¬ 
placements, respectively, and p the density profile of the star. 
Likewise, the rotation kernel can be expressed as follows: 

k 2 + 7(7 + lV - 2^/7 - rf\ pr 2 

KnAr) = ----- i - 

[A 2 + lit +\)rj 1 ~ 2^77 - T? 2 ] pr 2 dr 

\(£-r]) 2 + (J 2 + {-2)if]pr 2 

= —!---. (3) 

f 0 R [(£ - rj) 2 + A 2 + Z ~ 2 )r] pr 2 dr 


2. The surface rotation rate is positive. When combined with 
the previous assumption, this implies that the entire rotation 
profile is positive, which avoids onion-type structures with 
counter-rotating shells. 

We return to Eq. 0} and do an integration by parts of the right 
hand side: 



(4) 


where we have made use of the fact that K n f is unimodular, can¬ 
celled out one of the terms, and used the index i as shorthand 
for ( 77 , 7). If the rotation profile has a discontinuity at 77 /, then the 
integration by parts can be carried out on the domains [ 0 , 77 /] and 
[rd,R] separately: 


Sj 

1 -Q 


Q(R) + [Q(r-)-Q(r d + )|£'w(r)dr 



(5) 


We note that to be consistent with our first assumption, it makes 
more sense if Q( LP > £ 2 ( 7 -^). Similar formulas can be obtained 
for multiple discontinuities. At this point, we introduce a first 
type of normalised rotational splitting: 


s 


l-cr 


( 6 ) 


A first, and rather trivial, inequality is immediately appar¬ 
ent from either Eqs. © or (0. Indeed, since (p < 0, then the 
last term(s) on the right-hand side is positive. Consequently, this 
leads to: 

s' > (7) 

Such an inequality is obvious: the left-hand side is a weighted 
measure of the internal rotation rate, and the right-hand side the 
surface rotation rate. Given the assumption f * 0 , the surface 
rotation rate is necessarily smaller than the internal rotation rate. 

At this point, we introduce a first type of integrated rotation 
kernel: 

Ii(r)= f K,(r')dr'. (8) 

Jo 


Bearing in mind that rotation kernels are only defined for non- 
radial modes (t > 1), it is straightforward to see that K n j(r) is 
positive for all r values. Furthermore, K n j is unimodular, i.e., 

$*K n Ar)dr = 1. 

2.1. Rotation profiles with a negative gradient 

2.1.1. Inequalities for the full domain 

At this point, we introduce two assumptions concerning the ro¬ 
tation profile: 

1. The gradient of the rotation profile is negative < 0). As 
mentioned in the introduction, there are stars whe re this is 
not the case (ISchou et all Il998l iKurtz et all l2014lf . but we 
expect this to be true in many cases, especially in sub-giants 
and giants. 


Given that K, is unimodular, we have I JR) = 1. Furthermore, 
since K, is positive for all r and only zero in isolated points, the 
function /, is strictly increasing. Given that 7,(0) = 0 by con¬ 
struction, I Jr) is strictly positive for r > 0. A visual inspection 
of such functions for a set of modes, such as is illustrated in 
Fig. El shows that these functions tend to “line up” rather than 
cross each other. In other words, if IJjq) < l/nd for a given 7 - 0 , 
where i and j represent two modes, then I Jr) < I Jr) for all r val¬ 
ues. A similar behaviour certainly does not apply to the kernels 
themselves, hence the reason why we work with the integrated 
kernels. Now, it turns out that the integrated kernels do cross 
frequently (see Fig.©, but the general trend still leads us in the 
right direction. Indeed, to make the argument more rigorous, one 
simply needs to find constants, a and h, such that the following 
inequalities hold: 

Vr e [0,/?], 


aljr) < I Jr) < hi Jr), 


(9) 
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The optimal values of a and b will then simply be: 


a = min 

re[0,R] 



b = 


max 

re[0,R] 


m ) 


( 10 ) 


Figure [2] shows the ratios of integrated kernels for two pairs of 
modes from which we determine a and b. Already, it is straight¬ 
forward to see that a < 1 < b by simply inserting r = R in 
Eq. 0. Furthermore, if a = 1 or b = 1, this implies that the in¬ 
tegrated kernels do not cross (although we do note that by con¬ 
struction, they will have the same values at r — 0 and r - R, i.e. 
0 and 1 , respectively). 



Fig. 1 . (colour online) Integrated normalised kernels of dipole 
modes for a 1 M e red-giant (model 1, see Sect. |3j. The lower 
curve corresponds to a more p-like mode where as the top curves 
are for more g-like modes. 


The behaviour in r — 0 is more complicated. Indeed, one 
needs to consider the limit of Ij(r)/Ij(r) in Eq. (ITOl) when r —> 0. 
Hence, it is useful to know the behaviour of the rotation kernels 
when r goes to 0. For ( > 1, the vertical and horizontal displace¬ 
ment behave as 0(r f 1 ). Therefore, the rotation kernels scale as 
r 2( when r goes to 0 as can be seen from Eq. (0. However, there 
is one exception to this rule. Indeed, when r goes to 0, the ver¬ 
tical and horizontal displacements satisfy the relation: £ ~ £tj. 
If we replace £ by Cjj in Eq. 0, then we obtain the following 
expression for the numerator: 

(2fi -e-\) T] 2 pr 2 = (2£ + \ ){{ - 1 )rf pr 1 . ( 11 ) 

This expression is zero when l — 1. Hence, in dipole modes, the 
lowest order terms cancel out and one needs to consider higher 
order terms. If we return to expression 0 and substitute i = 1, 
then the second part of the numerator drops out and we are left 
with (£ - rj) 2 pr 2 . The next order term for (£ - 77 ) behaves as 0(r 2 ) 
since only even powers of r intervene in the displacement func¬ 
tions of dipole modes. When squared and multiplied by pr 2 . this 
leads to an 0(r 6 ) behaviour, rather than the 0(r 2 ) behaviour ini¬ 
tially expected. Finally, if we return to the general case, the inte¬ 
grated kernels will simply behave as 0(r lt+1 ), except for dipole 
modes, for which 7,-(r) = Oir 1 ). Therefore, when considering 
the limit /,(r)// ; (r), the modes i and j need to have the same £ 
value or else one needs to have £ — 1 and and the other £ — 3. 
Otherwise, the limit will either be 0, thereby leading to a — 0, or 
infinite, thereby leading to b = 00 . 



Fig. 2. Ratios of integrated kernels for dipole modes. The top 
panel shows the ratio between the integrated kernels of a g-like 
mode and a p-like mode, whereas the bottom panel is for two g- 
like modes. As can be seen in the lower panel, the ratio crosses 
the value 1 , which means that the two integrated kernels have 
crossed. 

If we now multiply Eq. @ by — ^ and integrate over [0, R], 
we obtain the following inequalities: 

a [s'. - Q.(R)] < [s' - Q(/?)] < b [s' - Q(R)], (12) 

where we have made use of the assumption ^ < 0 , and of 
Eq. 0, which assumes that 11 is continuous. Bearing in mind 
that aQ.(R) < 11(7?) < /;!!(/?), one can simplify the terms propor¬ 
tional to the surface rotation rate, thereby leading us to our final 
set of inequalities: 

as) < s' < bs). (13) 

Although simpler, this last equation is slightly less restrictive 
than the previous form. Hence, one should use Eq. (fl2l > if the sur¬ 
face rotation rate is known. If the rotation profile, Q, is discon¬ 
tinuous at r r j, it is still possible to obtain the inequalities given in 
Eq. ([lA^, provided that H( r 7> > fl(rf), which is consistent with 
^ < 0 as pointed out earlier. As will be discussed in Sect. 12.1.51 
observational error bars also need to be taken into account when 
applying the above inequalities. 

We note, in passing, that the quantities R - r, also obey the 

rR 

above inequalities, where r, = J rKj(r)dr. Indeed, an integra¬ 
tion by part yields: 

J r*R r*R 

rKj{r)Ar = R— I 7,-(r)dr. (14) 

0 Jo 

Hence, integrating Eq. 0 over [0, R] leads to: 

a(R - rj) < R - r, < b(R - rj). (15) 

If a or b is equal to 1 for a given pair of modes (i.e. if their in¬ 
tegrated kernels /, and I s do not cross, as noted above), then the 
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slope of line connecting (R - r,, s') to (R - rj, s') must be posi¬ 
tive. However, it turns out that integrated kernels do cross fairly 
frequently (as shown in Fig. [3}. so one cannot rely on the slope 
to decide whether a particular pair of modes obey the above in¬ 
equalities. Instead, one needs to apply the inequalities system¬ 
atically. We nonetheless expect there to be a general trend in an 
(R - rj, s-) diagram, based on these considerations. 



Fig. 3. Cumulative number of crossings between the 7,(r) inte¬ 
grated kernels as a function of position, r/R. These crossings 
were calculated using the £ — 1 modes from Model 1, which is 
introduced in Sect. 13.21 


Bearing in mind that the original goal, as described in the 
Introduction, is to be able to find rotation profiles which are de¬ 
creasing and which do not change signs, it is interesting to look 
at the above problem the other way around and see whether it is 
possible to construct such rotation profiles for given rotational 
splittings, assuming they obey the above inequalities. We begin 
by denoting y = 4. The quantity y is between a and b since the 

J 

rotational splittings satisfy the above inequalities. Furthermore, 
the function jj^- ranges from a to b, by definition of these con¬ 
stants. If this function is continuous, then there exists a point r y 
such that T- 2 - = y. It is then straightforward to see that the fol- 

1 jVy) 

lowing rotation profile reproduces the rotational splittings and 
satisfies the above constraints: 




if 0 < r < r y 
0 if r y < r < R 


(16) 


If the function never takes on the value y due to a discontinu- 
ity (for instance in a model with a discontinuous density profile), 
then one could define a rotation profile which is discontinuous 
at the points r\ and ro, defined such that < y < and 
solve the relevant system of equations to find by what amount 
the rotation profile changes at each discontinuity. A more seri¬ 
ous difficulty occurs if r y — 0 and no alternative r values could 
be used to construct a solution. In such a situation, one can only 
get arbitrarily close to the solution by setting a discontinuity at 
e, finding the corresponding rotation profile, then decreasing e. 
Hence, in summary, the above inequalities provide a necessary 
and nearly sufficient condition on the rotational splittings of two 
modes for the existence of a decreasing and positive rotation pro¬ 
file. 


Had the above inequalities not provided a nearly sufficient 
condition for the existence of such profiles, then one would be 
left wondering if a more stringent set of criteria could be de¬ 
duced. However, the above results suggest that these are the most 
restrictive conditions one could find for a given pair of modes. 
As we will, however, see in a later section, they do not provide 
the most restrictive conditions for a set of 3 or more modes, 
as they apply to 2 modes at a time. Also, one must not forget 
that even if rotational splittings do obey the above inequalities, 
it does not guarantee that the true rotation profile is indeed de¬ 
creasing everywhere. Indeed, for any set of rotational splittings, 
it is always possible to find rotation profiles which have a posi¬ 
tive gradient somewhere in the star and/or a sign change. 


2.1.2. Inequalities where the centre is excluded 

If an upper bound on - ^ is known in the most central regions 
of the star, then it is possible to derive slightly different inequal¬ 
ities which exclude these regions and potentially lead to tighter 
constraints. We start by defining new constants, a* and b*, as 
follows: 


a 


★ 


min 

re[r 0 ,fl] 


IM \ 


b* 


max 

re[r 0 ,R] 


IM \ 


(17) 


where ro is the upper bound of the central region under consid¬ 
eration. Hence, 


Vr e [r 0 ,7?], a*I dr) < 7,(r) < b*L(r). 


(18) 


We then multiply this equation by and integrate over [ro, 7?]. 
Remembering that 


C R dQ , p dQ 

1. dT ,,(rMr ‘ s ‘- n, R > + l 17 


Mr) dr. (19) 


where k corresponds to i or j, we finally obtain, after some rear¬ 
rangement and cancelling out the surface rotation terms: 


X ro flQ 

— \a*Ij(r) - 7,(r)] dr < s' 

S ' J + l ^[ b * Ij(r) ~ Ii(r) ] dr - 


< b*s' f 


( 20 ) 


At this point, we introduce B, the upper bound on — ^ over the 
interval [0, ro]. The first integral in the above expression can be 
bounded as follows: 


-r 

Jo 


f [a*7/r)-7,(r)]dr 

r° df! r * 1 

< - J — max 10, a Ij(r) - 7,(r)J dr 


< B 


fo fo, a*Ij(r) — 7,(r)] dr, (21) 


where we have also made use of the assumption < 0 over 
the interval [0, ro]. We note that, if we substitute a for a*, 
then max |(), alj(r) - 7,(r)j is always zero, thereby cancelling out 
the right-hand side. The integral term would then drop out of 
Eq. (l20l >. A similar manipulation with the second integral term 
leads to: 

£° f |^*//(r) - 7;(r) | dr 

X r ° r 

max |o, 7;(r) - 7>*7 7 (r)J dr. 


(22) 
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Once more, if we substitute b for b*, the right-had side cancels 
out. When substituted into Eq. (l20b . these inequalities lead to: 


- ® J o max \o,a*Ij(r) - /,-(r)] dr < s' 

n> o 

< b*s'j+!B I max [0,//(r)-i*/Xr)]dr. (23) 


The advantage of this equation over Eq. (Il3l) is that the constants 
a * and b* may be much more constraining than a and b, given 
that the latter may be overly affected by the very inner regions. 
In particular, this would allow comparisons between splittings 
for different t values, for which the constants a or b could be 
0 or infinite, respectively. The drawback is trying to find an ap¬ 
propriate value for S. One can expect '£ to become small in the 
central regions since ^ = O(r). 

In much the same way as was done above, one can check to 
see if the above inequalities provide a sufficient condition for ob¬ 
taining a positive, decreasing rotation profile, subject to the sup¬ 


plementary condition - ^ < £ over the interval 0 < r < iy. If 42 
is prescribed over the interval [0, r (l [, one could apply a similar 
approach to what was done above and attempt to define a radial 


coordinate r y * such that such that y* 


t+fo 0 fb(>)dr _ /,( V ) 

W)’ 

If such a point exists, then it is possible to construct a rotation 
profile which is decreasing and satisfies the rotational splittings. 
However, Eq. (1231) does not guarantee that y* is between a* and 
b*, and hence that /y exists. Only Eq. ( 1201 ) provides such a guar¬ 


antee. One may try to adjust the function ^ over the interval 
[0, ro[, but would only succeed in enforcing either a* < y* or 
y* < h*. One could then take on a different approach, and define 
the points r a + and tv such that yjfyyj = a* and j j = b*. As 


opposed to r y *, these points are guaranteed to exist. However, as 
described in Appendix [A] if one defines a rotation profile with 
discontinuities at those points, and which reproduces the rota¬ 
tional splittings, only Eq. (l20l > guarantees that the rotation profile 
will decrease across these discontinuities. Hence, only Eq. (l20l > 
provides a necessary and sufficient condition for the rotation 
splittings to correspond to a positive decreasing rotation profile, 
provided it is decreasing over the interval [0, rot- Equation (l23l >. 
in contrast, constitutes a necessary but insufficient condition on 
the rotational splittings. 


2.1.3. Three or more modes 

So far, we have only considered two modes in isolation. 
However, in typical stars, a larger number of rotational splittings 
are observed. Of course, one can always apply the above inequal¬ 
ities to every pair of rotational splittings. However, it is obvious 
that a more complicated strategy than what was described in the 
previous section is needed in order to construct a rotation pro¬ 
file which satisfies all of the rotational splittings simultaneously, 
and yet still has a negative gradient throughout the star. Indeed, 
certain properties of the rotation profile only emerge when a suf¬ 
ficient number of rotational splittings are used together. Hence, 
in what follows, we briefly explore a way of generalising the 
above inequalities to more than two modes. Let us consider an 
inequality of the following form: 

Vr 6 [0,/?], £ ai l,(r) < £ b jlj( r h (24) 

i j 

where we are assuming that the a, and bj are positive. Then us¬ 
ing the same methodology as above, we deduce the following 


inequalities on the corresponding rotational splittings: 

Y aiS 'i ~ Y bjS 'j- (25) 

' j 

We note that it was possible to remove the surface rotation rate 
because < Yijbj (as deduced from Eq. (l24l > for r = R). 
If one excludes the central region, [0, ro], then the following in¬ 
equality is obtained: 


nro 

I max 0, ^ a,/;(r) - ^ bjlfr) dr<^b 
; JO 


( 26 ) 

where £ is an upper bound on over the interval [0, ro]. One 
can also obtain the following condition by making use of the 
inequalities ^/ ; (r)dr < 0, thereby allowing us to remove 
these terms altogether: 


/,(r)drj < ^ b,s'. (27) 

Although less restrictive than Eq. (l26l >. this latter inequality has 
the advantage of being linear with respect to the coefficients a, 
and bj. 

The main difficulty in the above inequalities is finding the 
constants a, and bj in such a way as to provide tight constraints. 
For instance, one could choose a set of bj values, and search 
for aj values which maximise the left hand side of Eqs. 4251) 
or (l26l > while respecting Eq. (l24l) . Equation (l24l) could be ap¬ 
plied at each mesh point, thereby providing a set of N linear in¬ 
equalities on the coefficients a,-, where N is the number of mesh 
points. Maximising the left-hand side of Eq. ( l25l > would then 
require methods from linear programming such as the simplex 
algorithm. A similar strategy can also be applied to Eq. 427). In 
contrast, maximising the left-hand side of Eq. 4261) would require 
a method from mathematical optimisation, due to its non-linear 
form. In the next section, we carry out a comparison between 
these inequalities and inverse methods, thereby providing fur¬ 
ther insights into how to choose the coefficients a,- and bj. 


2.1.4. Link with inverse methods 

In order to understand the link between the above inequali¬ 
ties and inverse methods, we start with Eq. 424) and subtract 
2, aJiir) from both sides: 

Vre[0,/?], 0 <Y J b j Ij(r)-Y j a i Ii(r)= r K(r). (28) 

j i 

The right-hand side of this inequality is a linear combination of 
kernels that plays exactly the same role as an averaging kernel 
in inverse methods. There is, however, one key differencqj 'K(r) 
needs to be positive everywhere, as indicated by the equation, 
in order for the inequality to carry through the integration onto 
the rotational splittings. In contrast, the averaging kernels from 
inverse methods are not obtained with such a goal in mind, and 
therefore frequently take on negative values. 

This condition can be relaxed if one excludes the centre, or 
some other region(s) in the star. Indeed, in such a situation, the 

1 A second difference between the averaging kernel and the function 
7C is that the former is normalised so as to yield a proper average of 
the rotation profile. This, however, is a minor issue for the inequalities 
presented above, since these are not affected by the normalisation of 7C. 
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combined kernel, 'K, only needs to be positive in those regions 
which have not been excluded. The upper bound on -jp can 
then be used to constrain the supplementary terms which arise 
from the excluded regions, as is done, for instance, in Eq. (l26l >. 
Hence, one could, in principle, take a set of inversion coeffi¬ 
cients, exclude the regions where the averaging kernel is neg¬ 
ative, and deduce inequalities on the rotational splittings for a 
given upper bound, B. Conversely, one could search for the lim¬ 
iting value of B beyond which the inequalities break down. If 
the averaging kernel is well localised, one will have primarily 
tested whether the rotation gradient is negative in that particular 
region. This approach could also be combined with a simplex 
method, as described above, in order to find tighter constraints 
on the splittings. 

2.1.5. Error bars 

In practise, it will be necessary to take into account error bars on 
observed rotational splittings. As will be described in this sec¬ 
tion, this introduces complications when interpreting the above 
inequalities. Indeed, if one of the inequalities is not satisfied, one 
can only deduce that the rotation profile has a positive gradient 
with some probability that needs to be determined. In order to 
illustrate this, we start with Eq. (l25l > as a generic form for the 
inequalities, group the non-zero terms together on the left-hand 
side, and introduce error terms: 

Y aiS '> ~ Y bjS 'i + Y aie '‘ - Y bj£ 'j ~ °’ (29) 

i j i j 

s' e' 

where e' is the error realisation on a given normalised splitting, 
s' represents the combined splittings, and s' the combined errors. 
For a true violation of the inequality, one needs s' > 0. However, 
since only the measurement s' + s' is available, one has to eval¬ 
uate the probability that s' + s' > s'. An obvious approach is to 
compare the lcr error bar on s', which we will denote cr' , with 
the measurement. This can be obtained as a quadratic sum of the 
individual lcr error bars: 

cr' = af(cr') 2 + ^ bj(o-'j) 2 , (30) 

where the cr' are the lcr error bars on the individual normalised 
splittings. Hence, if s' + s' is equal to 3cr', and s' follows a nor¬ 
mal distribution, the probability that the inequality is violated is 
99.87 % (where we’ve taken into account the fact that there is vi¬ 
olation in only one of wings of the distribution for s'). It would 
be tempting to conclude that this is the probability that the rota¬ 
tion profile has a positive gradient. However, one will typically 
test a large number of inequalities which increases the chances 
of finding large deviations on one of the s' values and hence 
of having a false alarm. Given the correlations between the dif¬ 
ferent inequalities, it is not straightforward how to calculate the 
probability of a false alarm. As will be described in Sect. [3 we 
will carry out Monte-Carlo simulations to estimate such a prob¬ 
ability. 

2.2. Rotation profiles subject to Rayleigh’s stability criterion 


energy by redistributing its angular momentum ( e.g. IRieutordl 
ll997f) . Mathematically, this is expressed by the condition: 


d(m 4 Qr) 

dm 


> 0, 


(31) 


where m is the distance to the rotation axis. If we furthermore 
assume that the rotation profile only depends on r, this criterion 
then becomes: 


d(r 4 Q 2 ) 

dr 


dln|Q| 

dlnr 


(32) 


Hence, this criterion gives the maximum rate at which a rotation 
profile can decrease before the fluid becomes unstable. We note 
that such a criterion is incompatible with a sign change in the 
rotation profile, except at discontinuities. However, in realistic 
stars, viscosity, even though it is small, would remove true dis¬ 
continuities, thereby suppressing sign changes. If a discontinuity 
were present in the rotation profile, it would be stable only if the 
absolute value of the rotation rate increases outward across the 
discontinuity, which is the opposite from what was considered in 
the previous section (if we ignore sign changes). Finally, in real 
stars, this criterion is only one of the terms in the more general 
Solberg-Hoiland criterion for convective stability (e.g. iMaedeii 
l2009h . Hence there could be situations where an unstable rota¬ 
tion profile is stabilised by, say, a p gradient. 

At this point, we will assume that the rotation profile does 
not change signs. Accordingly, we will see the consequences of 
the inequality d( ' d/ n> > 0 rather than those of Eq. 1121 ). In or¬ 
der to derive inequalities on the rotational splittings, we start 
from Eq. <JT]), do an integration by part, and divide both sides by 





= R 2 Q.(R) ■ 


r R d (r 2 Q) 

r MDdr' 

Jo (r ') 2 

Jo dr 

f R ^d r' 

L Jo (r ') 2 


dr. 


(33) 


If a discontinuity is present, for instance at r = rj, one would 
obtain the following formula: 


(l-Qj^d r 


r ^d r 

= R 2 Cl(R) + r 2 [fl(rj) - H(r+)] Jo 


I 

i: 


d(r 2 Q) 
dr 


f £2G d r' 

Jo (r') 2 U ' 


d(r 2 G) 

dr 


f* M2 dr , 

Jo (r') 2 U ' 

r MBl dr ’ 

Jo Cr"') 2 


f R M2 dr > 

Jo (r') 2 U ' 


f 

do 

dr 

dr. 


■ R km 


dr 


(34) 


Of course, the difference G(rJ) - A,) would have to be neg¬ 
ative, in keeping with Rayleigh’s criterion. We note, in passing, 
that the rotation kernels, K,, behave as r" in the centre, where 
n > 4. Hence the ratio, Kj(r)/r 2 always goes to zero in the cen¬ 
tre. From now on, we will use the following notation to designate 
a second type of normalised rotational splitting: 


2.2.1. Inequalities for the full domain 

A different set of inequalities can be obtained from Rayleigh’s 
stability criterion. According to this criterion, the angular mo¬ 
mentum must increase with the distance to the rotation axis: oth¬ 
erwise, the fluid will be dynamically unstable and can free up 




(35) 


A first inequality can be found straight away, in much the 
same way as was done above. The second term on the right-hand 
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side of Eq. (1331) is negative, given the assumption dl ' d) n> > 0. This 
leads to the following inequality: 

St < R 2 Q(R). (36) 

This inequality sets a lower limit on the surface rotation rate 
based on the internal rotation rate, as measured by the rotational 
splitting. This is consistent with Rayleigh’s criterion which gives 
the maximum rate at which the rotation profile can decrease. 

We then introduce a second type of integrated kernel: 



Examples of this type of integrated kernel are shown in Fig. [4] 



Fig. 4. (colour online) Second type of integrated normalised ker¬ 
nels for the same modes as in Fig. [T] (the same colour scheme 
is used in both figures). These integrated kernels are very sim¬ 
ilar, so it is difficult to distinguish them, even though there’s a 
mixture of p-like and g-like modes. 

As was done above, we find, for a given pair of modes (i, j), 
inequalities of the following form: 

Vr 6 [0,R], ajj(r) < 7,(r) < bJj(r), (38) 



Fig. 5. Ratios of second type of integrated kernels for the same 
pairs of modes as in Fig. [2] 

Hence, we rearrange this inequality so as to make positive terms 
appear, thereby leading us to our final form: 

bsj -(b- 1)Q (R)R 2 < Si < asj + (1 - a)Cl(R)R 2 , (41) 

where 6Sj is the error bar on the normalised splitting, Sj. Had 
we worked with a discontinuous profile, the same inequalities 
would be obtained, as long as the rotation rate increases across 
the discontinuities (which is the opposite to what we assumed in 
the previous section). Otherwise, supplementary terms for each 
discontinuity where the rotation rate decreases would need to 
be included, but again we emphasise that such discontinuities 
would not satisfy Rayleigh’s stability criterion. 

In analogy with what was done above, we introduce the 
quantity r defined as follows: 


where 


a = min 

re[0,R] 


JM \ 

wr 


b = max 

re[0,fi] 



(39) 


Figure 0 illustrates ratios of the second type of integrated ker¬ 
nels, from which we can determine a and b. Once more, we ob¬ 
tain the relations a < 1 < b by inserting r = 1 into Eq. (l38l >. 
In order to avoid having a - 0 or h - oo, the integrated kernels 
J, and Jj need to have the same behaviour in the centre. Hence, 
the £ values of modes i and j, should be the same or should be 1 
and 3. 

We then multiply Eq. (l38l > by - d( ' d) n) and integrate over the 
interval [0,R]. This yields, taking into account the sign change: 


a [s, - fi(R)R 2 ] > [si - Q(R)R 2 ] > b [s, - Q(R)R 2 ], (40) 

where we have made use of Eq. (l33l >. This time, given the in¬ 
verted inequalities, the surface terms cannot simply be cancelled 
out and must therefore be kept. Furthermore, because of the in¬ 
equality given in Eq. (l36l >. the terms in Eq. (l40t are negative. 


Cr^dr 

JO r- 




(42) 


Carrying out an integration by parts leads to the following equal¬ 
ity: 




Ji(r)dr. 


(43) 


Combining Eq. (l43t with the inequalities in Eq. (l38l) then yields: 

a(R - rj) < R- r, < b(R - rj). (44) 


Hence, the quantities R - r, obey the opposite inequalities to the 
normalised splittings S;. Accordingly, if a = 1 or b = 1 (i.e. if the 
integrated kernels J,(r) and Jj(r) do not cross), then the slope of 
the line connecting ( R-r ,, ?,) to (R - rj, sj) is negative. However, 
the integrated kernels do cross frequently, as illustrated in Fig. [6] 
so one should instead look at the general trend in an (R - r,, j?,-) 
diagram. Only a systematic application of the inequalities will 
yield conclusive results. 
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Fig. 6. Cumulative number of crossings of between the /,(r) in¬ 
tegrated kernels as a function of position, r/R. These crossings 
were calculated using the £ — 1 modes from Model 1, which is 
introduced in Sect. 13.21 


One can then try to see if the above conditions are sufficient 
for the existence of a rotation profile which satisfies Rayleigh’s 
criterion. This time, in keeping with the inequalities, we must 
assume the surface rotation rate is imposed. As was done previ- 
ously, we introduce the ratio y = R iQ {R y~. and its corresponding 


radial coordinate, ry, defined such that y ■ 


Jj(y) 


. This radial co- 


j. (r) 

ordinate is guaranteed to exist provided the function y/- 2 is con¬ 


tinuous. We then define a rotation profile which is discontinuous 
at r=y and where d( ' d) n) = 0 everywhere else: 




yr if 0 < r < rj 
% if ry < r < R ' 

r~ y 


(45) 


The constants p and A are obtained by imposing the correct sur¬ 
face rotation rate and rotational splittings. This leads to: 


A = 


R-Q(R) - 

Ji(r?) 


p = R-Cl(R) 


R-Q.(R) 




(46) 


It is straightforward to show that A is positive, using Eq. (136) . 
However, it is not clear whether p is positive or not, meaning 
the above rotation profile might be negative in the centre. Also, 
we were not able to show that \p\ < \p + d| which is needed to 
ensure that the absolute value of O increases across the discon¬ 
tinuity, and hence that Rayleigh’s criterion is verified. Instead, 
the above rotation profile only satisfies the criterion d( ' d| £i) > 0, 
which is slightly less restrictive than the condition we obtained 
when we assumed the rotation profile doesn’t change signs in 
addition to Rayleigh’s criterion. Accordingly, if p is positive, 
then \p\ < \p + T| and Rayleigh’s criterion is nearly satisfied ( i.e. 
d ° d/ ~~ 1 > 0). Hence, this simple strategy is unable to produce a 
fool-proof solution solely based on the above inequalities and 
seems to indicate that these do not provide a sufficient, or nearly 
sufficient, condition for the existence of a rotation profile which 
satisfies Rayleigh’s criterion. 


2.2.2. Inequalities where the centre is excluded 

If an upper bound on d( ' df n> in the central regions is known, then 
it is possible to derive some slightly different inequalities, much 


stellar rotational splittings 


like in the previous section. We begin by introducing the quanti¬ 
ties a* and b *: 


a 


★ 


min 

re[r 0 ,R] 


JM\ 

Jj(r))' 


b* 


max 

re[) 0 ,R] 



(47) 


Then, following a similar reasoning as in the previous section, 
we finally obtain: 

nr o 

b* sj - (. b * -1 )Q.(R)R 2 - S I max [o, i,(r) - b*Jj(r)\ dr 

< ^ 

r r ° 

< a*sj + (1 - a*)Cl(R)R 2 + SJ max [o, a*Jj(r) - J,(r) J df48) 

where & denotes the upper bound on d( d) Q) over the interval 

[0, r () ]. We do note that to first order, d behaves as 2rQ in 
the centre. 


3. Applications 

3.1. Direct application of the inequalities 

A natural and straightforward application of the above inequal¬ 
ities is to see whether the rotation profile of a given star has a 
positive gradient and/or fails to respect Rayleigh’s criterion. In 
what follows, we will therefore test these inequalities in artificial 
situations to see how effective they are. 

We start by using Model S (iChristensen-Dalsgaard et akl 
1 996) and impose the following rotation profile: 

fl(x) = A + B tanh ^ j - C tanh ^j, (49) 

where A = B = 1.25, C = 1.5, Xo = 0.220, x\ = 0.432, 
and d = 0.05. This rotation profile and its derivatives are dis¬ 
played in the upper panel of Fig. [7] As can be seen, it has a 
positive gradient in the first 33% of the star, and does not satisfy 
Rayleigh’s criterion over a narrow region around 0.45 R. We then 
considered a set of 44 modes with harmonic degrees £ = 1 — 3 
from which we calculated “observed” rotational splittings. We 
assumed these splittings follow a normal distribution, and their 
lcr error bars are 0.0067 /./Hz, which is somewhat smaller than 
what is achieved with the Kepler satellite. 

A systematic application of the inequalities (fl3T> and ( 1231 ) to 
all pairs of modes with degrees £ - 1 to 3 did not reveal a neg¬ 
ative gradient. However, using Eq. (125), we were able to place 
tighter constraints on the splittings. We systematically applied 
this equation, in which each splitting was isolated to one side, 
then the other, of the inequality, and a simplex algorithirQ used 
to optimise the coefficients of the remaining splittings. This ap¬ 
proach allowed us to deduce lower and upper bounds on each 
splitting, and to show that the rotational splitting of the £ — 3 
mode with the lowest radial order (n = 13) is approximately 
4.99cr beyond the bound predicted by Eq. (125) . i.e. we showed 
that s' lcr' — 4.99. A naive interpretation of this would lead to 
the conclusion that the probability of a false alarm is 3.0 x 10 7 . 
However, as discussed in Sect. 12.1.51 this does not take into ac¬ 
count the fact that we are testing multiple inequalities, which 
increases the probability of having a large deviation on one of 

2 The simplex algorithm was downloaded from 
http://algs4.c s.princeton.edu/6 5reductions/Simplex.j ava.html 
and is described in lSedgewick & Wavne ( 2011 ). 



























9 


D. R. Reese: Inequalities on stellar rotational splittings 



0.0 0.2 0.4 0.6 0.8 1.0 

Position, r/R 

Fig. 7. (Colour online) (Upper panel) Rotation profile and as¬ 
sociated derivatives used in first test case on Model S. (Middle 
panel) Combined kernel from simplex method which led to a 
detection of a positive gradient. (Bottom panel) SOLA inversion 
results on ^ along with lcr error bars. 


them. We therefore carried out a Monte Carlo test with 10 6 real¬ 
isations of the splittings. We calculated the maximal value of 
s' /cr' for each realisation, using the coefficients from the in¬ 
equalities based on the mode pairs and based on the simplex ap¬ 
proach, and counted the number of times the threshold 4.99 was 
exceeded. This occurred in 166 cases, thereby implying a proba¬ 
bility of 1.7 x 1(T 3 4 for a false alarm, i.e. substantially larger than 
the above value. A similar test, using splittings from a decreas¬ 
ing profile, yielded a probability of 1.5 x 10 4 , thereby showing 
the result is robust. Even then, this new value is likely to be an 
underestimate, because the Monte Carlo simulation was carried 
out for fixed coefficients. Ideally, one would want to repeat the 
simplex calculation for each realisation of the splittings, since 
the results are optimised for the given set of splittings. However, 
the numerical cost of such an approach is too high to obtain a 
result in a reasonable amount of time. 

The reason why the inequality failed for the £ = 3 mode is 
straightforward to understand. The middle panel of Fig. |7] shows 
the resultant combined kernel, 7C. This kernel has a maximum 
amplitude where the rotation profile increases sharply and a low 
amplitude near the sharp decrease. Hence, the integral of the 
product K is negative which leads to a break-down of the 



2500 3000 3500 2120 2160 

Frequency (in p, Hz) 


Fig. 8. (Colour online) Rotational splittings and lcr error bars 
(solid red lines), constraints deduced from inequalities on pairs 
of modes (dotted blue lines), and constraints from multiple 
modes using simplex method (dashed green lines). The lighter 
lines denote the lcr error bars. Right panel : zoom oil £ — 3 mode 
where an inequality breaks down by 4.99cr (as indicated by the 
arrow). 


inequalities for sufficiently small error bars. The lower panel of 
Fig-Hshows the results obtained from a SOLA inversion of 
Given that only the £ — 1 to 3 modes were used, the averaging 
kernels are poorly localised and the resultant inverted profile a 
poor match to which inspires little confidence in the inver¬ 
sion. Nonetheless, the inversion does detect a positive gradient 
with a 2.87cr margin. Hence, the above inequalities, in conjunc¬ 
tion with the simplex method, have yielded firmer results than 
the inversion. 

In contrast to the tests based on finding a positive gradient, 
those based on Rayleigh’s criterion failed to detect a decrease in 
the angular momentum. Possible explanations for this failure is 
that the region with this decrease is either too narrow, or poorly 
located, given the set of modes we considered. 

We also considered a second rotation profile with model S: 


Q(x) = 


A B 


(50) 


where x = r/R, A = 1 pHz, B = 0.5/tHz, C = 0.3, D = 0.4, 
and E = 0.2. This rotation profile has been designed to have a 
positive gradient in a localised region, fail Rayleigh’s criterion 
in another region, and have a zero derivative in the centre, as 
expected from regularity conditions. Figure |9]illustrates this ro¬ 
tation profile, as well as df2/d.r and d(x 2 Q)/dx. 

A set of 282 modes with £ ranging from 1 to 20 were calcu¬ 
lated using the ADIPLS pulsation code (IChristensen-Dalsgaardl 
120081 ). Such an extensive set is naturally out of reach for stars 
other than the sun, given cancellation effects, but is what is 
needed to carry out a sufficiently detailed inversion of the ro¬ 
tation profile. Rotational splittings were calculated using the 
above profile, and their lcr bars were set to 0.0033 pHz, which is 
smaller than what is currently achieved with the Kepler mission. 

Figure m shows the result of a SOLA inversion 
(Piiners & Thompson! 1199^ . 1994 ) of the gradient of the rota¬ 
tion profilfl The lcr error bars are displayed in this plot. The 
positive gradient shows up clearly. 


3 In order to cancel the surface term, we included the supplementary 

constraint K m$ (R) = - 2, c,/,(/?) = 0 using a Lagrange multiplier. This 
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Radius, r/R 

Fig. 9. (Colour online) The rotation profile, as based on Eq. (l50l >. 
and its derivatives. 



Position, r/R 

Fig. 10. (Colour online) Inverted rotation profile gradient. The 
continuous line shows the actual rotation profile gradient, 
whereas the blue crosses indicate inversion results and associ¬ 
ated lcr error bars. 

In contrast to the inversions, a systematic application of the 
inequalities (IT3l > and (l23l) to all pairs of modes did not reveal a 
positive gradient. The same also applied to tests based on mul¬ 
tiple modes, as described above (we limited ourselves to modes 
with £ — 1 to 10, given the numerical cost of the simplex al¬ 
gorithm). We also applied a heuristic approach which consisted 
in recuperating the inversion coefficients and adjusting them so 
as to have a positive combined kernel everywhere. This strategy 
was to no avail even if the error bars were set to nearly 0. A 
possible explanation for these failures is the fact that the rotation 
profile decreases only over a very narrow region. 

One may also see whether it is possible to show that the un¬ 
derlying profile does not satisfy Rayleigh’s criterion. Figure HU 
shows inversiorQ results for the gradient of the angular momen¬ 
tum, using the integrated kernels defined in Eq. (ITTl) . As can be 
seen for the full set of modes (upper panel), the inversion is able 
to detect the region which does not satisfy Rayleigh’s criterion. 

We then systematically applied the inequalities (|4TT > for each 
pair of modes. Figure fl2l shows the maximal values of s' lcr' 
as a function of the harmonic degree. Three different cases are 
investigated. The solid line shows what happens if we use the 
exact surface rotation rate, Cl(R) = 0.0836 //Hz. Deviations big- 


can, however, lead to inversion coefficients which are not as optimal for 

the above inequalities, since these only require 'K(R) > 0. 

4 Once more, we removed the surface term by setting K mg (R) = 
- 2, CiJj{R) = 0 through a Lagrange multiplier. 



Position, r/R 

Fig. 11. (Colour online) Inverted profile for the gradient of the 
angular momentum, using £ — 1 to 20 modes (upper panel) and 
only £ = 8 modes (lower panel). The continuous line shows the 
actual gradient, whereas the blue crosses indicate inversion re¬ 
sults and associated lcr error bars. For the sake of clarity, the 
horizontal error bars (obtained from the averaging kernels) have 
been omitted in the lower panel. 

ger than 4<x start to appear at ( = 8. The false alarm probability, 
based on a Monte Carlo simulation is below 10 5 ( i.e. none of 
the 10 5 realisations produced a large enough value of s' lcr'). 
The dotted line corresponds to the case where the overestimate 
£l(R) = 0.5 //Hz is used. This time, the deviations exceeded 4<x 
for £ > 17. The false alarm probability is 1.4 x 10 '. Finally, the 
dashed line shows what happens when the centre is excluded, 
using ro = 0.01 and S/R = 1.0//Hz (which is nearly 50 times 
too large), and still keeping Q (R) = 0.5 //Hz. Deviations greater 
than 4cr appear at ( = 9. The false alarm probability turned out 
to be 1.3 x 10~ * 3 . We also carried out a Monte Carlo test with 
10 6 realisations in the first case (i.e. with the exact surface ro¬ 
tation rate, over the whole domain) using only £ = 8 modes. 
This yielded a false alarm probability of 2.1 x 10 4 (we note that 
the maximum value of e'/cr' is 4.61 in this case). In contrast, 
an inversion based on just these modes yields rather poor and 
inconclusive results due to poorly localised averaging kernels as 
illustrated in the lower panel of Fig. [IT](we do note that one of 
the points in the inversion does suggest a decrease in the angu¬ 
lar momentum, but even the lcr error bar allows for a positive 
gradient). 

One can also investigate the use of these inequalities in red 
giant stars. However, these seem to be less effective for such 
stars, since the rotation kernels are mainly sensitive to the core 
or surface regions and hardly in between. Accordingly, it is not 
too difficult to construct pairs of rotation profiles, one of which is 
decreasing the other not, and which produce essentially the same 
rotational splittings (i.e. to within some error bar). Figure [13] 
shows such a pair - the corresponding rotational splittings, when 
calculated in Model 1 (or Model 2 - these models are intro- 
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Harmonic degree, i 


Fig. 12. Maximal values of s'/cr' as a function of f: for inequali¬ 
ties based on Rayleigh’s criterion. Three different test cases are 
investigated (see text for details). 


lows, we will show a situation where a slight mismatch between 
an “observed” star and a reference model leads to one of the 
inequalities breaking down, in spite of having a rotation profile 
where the gradient is negative everywhere. To do so, we start off 
with two very close 1 M 0 models of red-giant stars. The second 
model is the same as the first, except for an ad-hoc 50% decrease 
of the density near the surface. This is achieved by multiplying 
the density profile by the following function: 


A{l+tanh[^=^ll 

f(r) = 1 + —- „ 1 Jf 


(51) 


where rj = 0.995, A — 2 x 10 3 , and A = -0.5. Table Q] 
gives the characteristics of these two models. We note in pass¬ 
ing that the mass difference between these two models corre¬ 
sponds to 73% of Mercury’s mass. Figure fl4l shows the inertia 
of their t - 0 and 1 modes as a function of frequency, where 

In,( = / 0 M ||£|| 2 d m/[MZ 2 (R)l 


duced in the following section), are within 0.01 pHz. Indeed, 
as explained in iGoupil et al.l d2013t) , the rotational splittings of 
red-giants are primarily sensitive to the average core and surface 
rotation rates. Obviously, the above inequalities will not detect 
a positive gradient in such a situation, since there exists at least 
one decreasing rotation profile which reproduces the splittings 
within the error bars. Including higher ( modes doesn’t seem 
like a viable solution, as these are harder to detect and will expe¬ 
rience a stronger trapping, either in the p- or g-mode cavities. In 
order to cause one of the inequalities to break down, one needs 
to have a rotation profile which is increasing throughout most of 
the star. However, this will probably cause the splittings of the 
p-like modes to be comparable or larger than those of the g-like 
modes - already a very strong indication of an increasing rota¬ 
tion rate prior to any tests or inversions. In the following section, 
we nonetheless show how the above inequalities can be useful 
when studying red giants. 



Table 1. Characteristics of the two models. The second column 
gives the characteristics of Model 1, whereas the third column 
gives the relative difference between Models 2 and 1. 


Quantity 

Model 1 

Relative difference 

Mass (in M 0 ) 

1.000 

-1.22 x 10-' 

Radius (in R Q ) 

4.107 

0 

y/GM/R 3 (in pHz) 

12.00 

-6.1 x 10- 8 



Frequency, v (in /j, Hz) 


Fig. 13. (Colour online) Rotation profiles which give the simi¬ 
lar splittings (to within 0.01 pHz) in Model 1 (from Sect. 13.21) . 
Profile 1 is also used to calculate the “observed” rotational split¬ 
tings (using Model 2) in Sect. 13.21 


3.2. The effects of a mismatch between an observed star 
and the reference model 

At this point a natural question arises: if the above inequalities 
break down is this necessarily caused by a positive gradient in 
the rotation profile, or can there be other causes? In what fol- 


Fig. 14. (colour online) Inertias of { = 0 and 1 modes in models 1 
and 2 as a function of frequency. 

In what follows, the frequencies of Model 2 will be treated as 
the observations, and Model 1 will be used as a reference model. 
Furthermore, the “observed” frequency spectrum will include all 
of the ( — 0 modes over a given frequency range, and only a 
subset of the t — 1 modes. Indeed, only the most p-like dipole 
modes are expected to be visible in observed stars. We therefore 
selected the dipole modes at local minima of the inertia curve 
as well as the adjacent 2 modes on either side. Figure fl5l shows 
an echelle diagram with the frequencies from both models. We 
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assume lcr error bars of 0.01 //Hz on the frequencies, which is 
typical for the Kepler mission. Table [2] gives the value ofyj? ed = 

1 M l y ref —y obs \” 

Ji E/=i I - —' I where M is the number of frequencies. This 

quantity characterises the departure of the frequency spectrum 
of Model 1 from that of Model 2. 


180 

„ 160 
N 

K 

^ 140 

o 120 

0) 

l=j 

& 100 

QJ 

(l 

fc. 

80 


0 5 10 15 

Reduced frequency (in /./,Hz) 

Fig. 15. Echelle diagram comparing the frequencies of Model 1 
(with and without the KBCD surface corrections) and Model 2. 
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Table 2. Seismic y^ ed values for Model 1, with and without the 
KBCD surface corrections. 


Frequency set 

Wed 

Model 1 

1364.5 

Model 1 + KBCD 

19.4 


As can be seen from Table [2] the pulsation spectrum of 
Model 1 is a poor fit to the observations. Our observational error 
bars would therefore exclude Model 1 immediately. However, 
surface corrections based on iKieldsen et al.l (20081) (KBCD, 
hereafter) are frequently included in seismic studies of observed 
stars. We therefore applied such corrections to the frequencies of 
Model 1, using the following formula: 


corrected 

y n,e 


Vn,t 


( £ \^KBCD 

fl KBCD(— j 

Qn.e 


(52) 


where v n / are the original frequencies, y“j Tected the corrected fre¬ 
quencies, vo a reference frequency, «kbcd a multiplicative con¬ 
stant, /; K bcd the exponent used in KBCD’s power law, and Q n j 
the mode’s inertia, divided by the inertia of an £_= 0 mode inter¬ 
polated to the same frequency (e.g. lAerts et al.L[201OL page 470). 
We choose the values vo = 125 //Hz, which is in the middle of 
the frequency range, c/kbcd = 0.40315 in order to reproduce 


as accurately as possible the frequency deviations in Model 2, 
and (?kbcd = 4.9, the value obtained in KBCD after calibrat¬ 
ing Model S to the sun. The division by Q n j is a later addi¬ 
tion to the formula, which proved to be necessary in the case of 
red giants, due to the large variation s of the inertia of non-radial 
modes ( e. g. iDeheuvels et all 120 1 2l) . The second line of Table [2] 
shows the impact of including this surface correction. Although 
the isstill above 1 and is worse than the best fitting models in 
Deheuvels et al. ( 20 T2h, i t is better than any of the best fits found 


Deheuvels et al. ( 1201 41. Hence, such a model would probably 


be retained as a viable reference model with which to study the 
“observed” star (i.e. Model 2). 

We calculate rotational splittings in Model 2, using the fol¬ 
lowing rotation profile: 


O (r) — A + 


B 

1 + C(r/R) 2 


(53) 


where A = 0.1 //Hz, B = 0.9 //Hz and C = 9 x 10 3 . We note that 
this rotation profile decreases with r, but still satisfies Rayleigh’s 
stability criterion. This rotation profile is illustrated in Fig. [13] 
(see Profile 1). The resultant rotational splittings for model 2 are 
shown in Fig. [16] and these will serve as our observations. We 
assume lcr error bars of 0.007 //Hz on these splittings, which 
is approximately V2 times smaller than the error bars on the 
frequencies, which is what you would expect from generalised 
£ = 1 splittings. 



Fig. 16. Rotational splittings in model 2. The diamonds corre¬ 
spond to the splittings which were used in this study, whereas 
the crosses represent those which were discarded. 


Then, we systematically apply Eq. (IT3l) to every pair of 
dipole modes. For each mode, we end up with N - 1 upper and 
lower bounds, where N is the number of rotational splittings. 
Figure[l7](upper panel) shows the bounds which are obtained for 
one of the p-like modes in Model 1. According to this approach, 
this splitting is compatible with the lower and upper bounds from 
the other modes to within the 1.22cr. However, if the centre is ex¬ 
cluded, as explained in Sect. 12.1.21 then it is possible to find an 
inequality which breaks down by 5.87<x. The middle panel of 
Figure [17] shows the bounds obtained for the same mode, but 
using Eq. [23] with r 0 = 2 x 10~ 4 R and RB - 5 x 10 5 //Hz. 
We note that the actual upper bound over the interval [0, r (l ] is 
R'B = 3.270 //Hz, i.e. more than 5 orders of magnitude smaller. A 
Monte-Carlo simulation based on the coefficients deduced from 
inequalities on mode pairs found three cases where e'!<j' > 5.87 
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in 10 6 realisations, thereby indicating a false alarm probability of 
3 x 10 6 . The lower panel shows what happens when applying 
the same tests using Model 2 as a reference model. This time we 
set RB = 3.270pHz to be as restrictive as possible. Even then, 
none of the inequalities break down, which is consistent with the 
fact that we are using the true model and the rotation profile is 
decreasing. 
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Fig. 17. (colour online) Upper and lower bounds on a given split¬ 
ting, deduced from the N— 1 other rotational splittings, assuming 
that the rotation profile decreases. In the upper and middle pan¬ 
els, Model 1 has been used as the reference model, whereas in 
the lower panel. Model 2 is the reference model. In the lower 
two panels, a more restrictive set of inequalities is obtained by 
excluding the centre, as described in Sect. 12.1.21 The frequency 
of the original mode is shown by the red vertical dotted line, 
whereas the blue horizontal dotted and solid lines indicate its ro¬ 
tational splitting plus or minus the error bar. The dark blue solid 
vertical segments indicate the upper and lower bounds deduced 
from the other modes, and the lcr error bars are indicated by the 
light blue extensions. 
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It is relatively straightforward to understand why Model 1 
leads to a break down of one of the inequalities. Indeed, as can 
be seen in Fig. [141 the mode inertias of Models 1 and 2 do not 
match very well, especially for the high frequency p-like dipole 
modes. This modification of mode inertias stems from the fact 


that surface effects act differently on p and g-modes and will 
therefore shift the former with respect to the latter. This, in turn, 
modifies the avoided crossings which occur between them, and 
hence, the inertias of the resultant modes. Such an effect is likely 
to be strongest in red giant stars where the spacing between 
consecutive g-modes can be comparable to the frequency shifts 
caused by surface effects. In contrast, applying a surface cor¬ 
rection recipe like the one described in KBCD merely corrects 
the frequencies and not the inertias. This leads to inconsistencies 
between the mode inertias and period spacings, as illustrated in 
Fig. m an d could mask mismatches between observations and 
what seems to be the best-fitting model(s). 



50 60 70 80 90 100 

Period spacing (in s) 


Fig. 18. Normalised mode inertias as a function of averaged pe¬ 
riods spacings (defined as 7i„ = (l/v„_i - l/v„ + i)/2, where n is 
the radial order). The lower panel corresponds to Model 1 with 
a 90 % decrease in density near the surface, whereas the upper 
panel uses Model 1 with KBCD surface corrections, adjusted so 
as to match (as best as possible) the frequencies of the modified 
model. The increased scatter at low inertias in the upper panel 
illustrates the mismatch between period spacings and mode in¬ 
ertias. 


One can also apply the tests based on Rayleigh’s criterion. 
We therefore systematically applied Eq. ( ITO using the exact sur¬ 
face rotation rate, to every pair of modes. However, this set of 
inequalities turns out to be far less stringent than the first set of 
inequalities. Accordingly, none of them break down, even when 
the centre is excluded, as described in Sect. 12.2.21 using the most 
restrictive value of B possible. FigurelT9lshows lower and upper 
bounds for the rotational splitting of the dipole mode (n, l) = 
(-61,1) with r 0 /R = 2x 1(T 4 and B/R = 4.037 x 10~ 4 pHz. 

One could then turn the inequalities based on Rayleigh’s cri¬ 
terion the other way around and search for a lower bound on 
the surface rotation rate. Equation (l36l > already provides a first 
condition. Applying this condition using Model 1 as the refer¬ 
ence model leads to Q (R) > 6.54 x 10~ 7 pHz whereas Model 2 
(the true model) leads to Q.{R) > 6.04 x 10 7 pHz. The true sur¬ 
face rotation rate is more than 5 orders of magnitude larger, with 
Cl(R) = 0.1001 Op Hz. However one can also rearrange Eq. (l48l > 
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Fig. 19. (colour online) Same as Fig. [17]but using the inequali¬ 
ties based on Rayleigh’s criterion. The centre has also been ex¬ 
cluded, as described in the text, thereby leading to more stringent 
conditions. The horizontal red dotted line near the top of the fig¬ 
ure is an upper limit deduced from the surface rotation rate. The 
error bars which extend off the bottom of the plot go to negative 
values. 


to obtain a more restrictive condition: 


b*Sj - Si - sfg max [o, J,(r) - b*Jj(r)\ dr 
R 2 (b * - 1) 

si - a*Sj - Sf 0 '°max [o, a*Jj(r) - J t (r) J dr 
R 2 (l - a*) 


< fi(R) 

< Cl(R) 


(54) 

(55) 


Applying these conditions using ro = 2 x 10 4 and 'BjR — 
4.037 x 10 -4 /rHz yields £l(R) > 4.89 x KTVHz for Model 1 
and £l(R) > 1.52 x 10~ 5 yuHz for Model 2, which is still 3 to 4 
orders of magnitude below the true value, and furthermore, very 
model dependant. Hence, even this way around, the inequalities 
based on Rayleigh’s criterion only provide very weak constraints 
and one would probably be better off either trying to measure the 
rotation rate through photometric modulation induced by spots 
(if present) or through spectroscopy. Alternatively, one could ap¬ 
ply the relation derived bv iGounil et al.1 ( 2013 1 which relates the 
ratio between the average surface and core rotation rates of red 
giants to the slope of the rotational splittings as a function of if, 
the ratio of the kinetic energy of the mode in the g-mode cav¬ 
ity to the total kinetic energy. Such a relation has already been 
tested on 6 subgiants observed by Kepler and shown to agree , 
in mo st cases, with results based on inversions dDeheuvels et all 
» . 


4. Discussion 

In this article, we have derived various inequalities which can 
be used to check whether a rotation profile has a negative gra¬ 
dient or satisfies Rayleigh’s stability criterion, for a given refer¬ 
ence model. We explored how these inequalities may be used in 
practical situations by investigating various test cases. In the first 
test case, a rotation profile which satisfies neither the assumption 
^ < 0 everywhere nor Rayleigh’s_criterion was imposed on 
Model S (iChristensen-Dalsgaard et all 119961) . A set of 44 split¬ 
tings for modes with C = 1 to 3 were used in the various tests. 


The above inequalities were able to detect the presence of a pos¬ 
itive gradient, but were unable to show a violation of Rayleigh’s 
criterion. Although an inversion also suggested the presence of 
a positive gradient, the resultant conclusions remained less firm 
than those derived by the inequalities. A second test case, also 
based upon Model S, used a different profile which also con¬ 
tained a region positive gradient on the rotation profile and one 
where the angular momentum is decreasing. This second case 
used a set of 282 modes (with £ — 1 to 20), thereby enabling a 
clear detection of both of these regions using inversions. In con¬ 
trast, the above inequalities failed to detect the positive gradient 
in the rotation profile, probably as a result of the narrowness of 
this region. However, these were able to detect a violation of 
Rayleigh’s criterion, even in situations where inversions yielded 
poor results, provided modes with £ > 8 were included. We ex¬ 
pect that this lower limit on the £ value is dependant on both the 
way in which the rotation profile violates Rayleigh’s criterion 
and on the size of the error bars on the rotational splittings. 

We also looked at other ways in which the inequalities may 
break down. For this purpose, we considered two models of red- 
giants which were identical except for an ad-hoc near-surface 
modification of the density profile in the second model. The first 
model was then used as a reference model in trying to inter¬ 
pret the rotational splittings of the second model, which were 
based on a rotation profile which has a negative gradient every¬ 
where and which obeys Rayleigh’s criterion. It turned out that 
one of the inequalities based on the rotation profile’s gradient 
was violated, whereas all of the inequalities based on Rayleigh’s 
criterion were satisfied. The reason for this is that the mode iner¬ 
tias differ in the two models, especially at high frequencies. This 
may come as a surprise since only the outer 0.5 % or so is dif¬ 
ferent between the two models. However, as argued above, near¬ 
surface effects shifts the frequencies of p-like modes relative to 
those of the g-like modes, thereby modifying the avoided cross¬ 
ings which occur between the two, and hence the resultant iner¬ 
tias. Frequency correction recipe^ for near-surface effects, such 
as the one presented in iKieldsen et al. ( 2008 ), are able to mask 
the frequency differences between the two models, but they do 
not correct the mode inertias, thereby leading to erroneous con¬ 
clusions on the rotation profile as shown in this test case. 

This leads to the natural question as to what to do if one 
or several inequalities are violated by a set of observed rota¬ 
tional splittings. In answer to this question, one needs to consider 
the following options. First, a different, more suitable, reference 
model needs to be found. We expect that in the case of red giants, 
the above inequalities could be quite constraining, due to very 
different mode trappings, and may therefore provide an indirect 
way of testing the structure of the reference model. Furthermore, 
such stars are the most likely candidates to having a rotation pro¬ 
file which decreases outward due to the spin-up of the core as it 
contracts and the slowing down of the envelope as it expands. Of 
course, a direct comparison of the non-radial mode frequencies 
will very likely provide tighter constraints, but so far, no one has 
yet fitted a red-giant model to within 1 cr of a set observed fre¬ 
quencies. Hence, if one adopts less stringent constraints on the 
selection of the model, and/or if frequency surface-corrections 
are included, the above criteria will be useful in determining be¬ 
forehand whether or not the chosen reference model is compat¬ 
ible with a decreasing rotation profile. The second option, is to 
consider that the rotation profile doesn’t have a negative gradient 
throughout the star and/or violates Rayleigh’s stability criterion 
(depending on which inequalities are not satisfied). As was al¬ 
ready stated in the introduction, there is convincing evidence of 
stars where the gradient of the rotation profile is positive, at least 









































D. R. Reese: Inequalities on stellar rotational splittings 


15 


in localised regions (ISchou et allll998llKurtz et alll2014t) . This 
raises the open questions as to how often such a situation occurs, 
and whether even a red giant couldn’t have a localised region 
with a positive rotation gradient, for example, in the outer con¬ 
vective envelope. However, one must not forget that the above 
criteria are not necessarily the most effective at detecting a ro¬ 
tation profile with a positive gradient. Accordingly, a mismatch 
between the reference model and the observed star may be a 
more likely explanation, which brings us back to the first op¬ 
tion. A last option would be to consider that the rotation profile 
is too rapid to be described by first order effects, thereby violat¬ 
ing the assumptions in this work. In such a situation, it would 
no longer be possible to apply linear inversion theory, except 
perhaps on generalised rotational splittings (defined as the fre¬ 
quency difference between prograde modes and their retrograde 
counterparts) o f acoustic modes in stars with nearl y uniform ro¬ 
tation profiles ( Rgesejstak, 20091). Recent wor ks dReese et al.L 
2006, iBallot et al. , 20101 Ouazzani et all 12013h have, however, 
shown up to what point a first order description of the effects of 
rotation are valid, and can therefore be used to decide whether 
or not the criteria described in this paper are applicable. 
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Appendix A: Attempting to construct a decreasing 
rotation profile when Eq. (23) is satisfied 

We start by introducing the following radial coordinates: 


= arg min 

re[r 0 .R] 



r b * = arg max 

re[r 0 ,R] 


IM \ 

IP-)}' 


(A.l) 


and define a profile as follows (where we have assumed r a * < 

n*y. 

{ A + n if r 0 < r < r a * 

Cl(r) = ifj. if r a * <r<r b + . (A.2) 

[0 if av < r < R 

A similar solution can be obtained if r b * < r a +. In this definition, 
we have deliberately postponed defining the rotation profile over 
the interval [0, ro[, except that we will impose its continuity at r {) 
(except if ro = r a *). Inserting this expression into Eq. ( IT9l > then 
yields the following system: 


, C" df2 

s, + J —Ij(r)dr = AIi(r a *) +/ili(r b *), (A.3) 

f'° dQ 

s'j + J —Ij(r)dr = Alj(r a *)+nlj(r b +), (A.4) 


the solution of which is: 


1 1 

f^-4+rf [b*w-iM] 


ijM | 

b* - a* 

] 

1 I 

+ [lM-a*Ij(r)} 

dr] 

iM | 

{ b *~ a * J 


(A.5) 

(A.6) 


From this we deduce that Eq. (l20b is a sufficient condition to 
ensure that A and /r are positive, and hence that £1 is decreasing. 
In contrast, Eq. ( l23l) . does not ensure that A and p are positive. 
To see this, we re-express the numerator of Eq. (lA.5b as follows: 

* , , f r ° dO r + n 

b Sj - Sj - J — max 10,7,(r) - b Ij(r) J dr 


~f min I 0 ’ dr ■ 

ii 

In order to ensure that the term I is positive according to Eq. (l23l >. 
one needs to set ^ to -£B over the interval where Ij(r)-b*Ij(r) is 
positive. Term II is necessarily smaller or equal to zero. In order 
to cancel it out, ^ needs to be zero over the interval where 7,-(r)— 
b*Ij(r) is negative. Likewise, in order to ensure that /./ is positive, 
one would need to make sure ^ = -£B when a*7,(r) - Ij(r) is 
positive, and ^ = 0 otherwise. However, there is no guarantee 
that these two sets of conditions are compatible. 


























